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ABSTRACT 

We present a sample of 485 photometrically identified Type la supernova candidates mined from the first three years of 
data of the CFHT SuperNova Legacy Survey (SNLS). The images were submitted to a deferred processing independent 
of the SNLS real-time detection pipeline. Light curves of all transient events were reconstructed in the qm, tm, *m 
and zm filters and submitted to automated sequential cuts in order to identify possible supernovae. Pure noise and 
long-term variable events were rejected by light curve shape criteria. Type la supernova identification relied on event 
characteristics fitted to their light curves assuming the events to be normal SNe la. The light curve fitter SALT2 was 
used for this purpose, assigning host galaxy photometric redshifts to the tested events. The selected sample of 485 
candidates is one magnitude deeper than that allowed by the SNLS spectroscopic identification. The contamination by 
supernovae of other types is estimated to be 4%. Testing Hubble diagram residuals with this enlarged sample allows us 
to measure the Malmquist bias due to spectroscopic selections directly. The result is fully consistent with the precise 
Monte Carlo based estimate used to correct SN la distance moduli in the SNLS 3-year cosmological analyses. This paper 
demonstrates the feasibility of a photometric selection of high redshift supernovae with known host galaxy redshifts, 
opening interesting prospects for cosmological analyses from future large photometric SN la surveys. 

Key words. Supernovae: general - Cosmology: observations 



1. Introduction 

The accelerated expansion of the Universe was revealed 
by studies of high redshift Type la supe rnovae (SNe 
la) (|Riess et aLl . 119981 : iPerlmutter et all . I1999D . Since then, 
second-generation experiments with efficient observation 
strategies were set up to achieve more precise cosmologica l 
measurements from SNe la, such as the SNLS (jAstier et aLl . 



[2006I) . ESSENCE (iMiknaitis et all. [20071) and SDSS-II 
Supernova Survey projects ( Frieman et aLl . [2"008[ ) . In order 
to meet the precision required by such measurements, sam- 
ples of SNe la rely on photometry to detect transient events 
and to measure their light curves, while SN la types and 
redshifts are provided by follow-up spectroscopy. In this 
paper, we describe a method based on photometry alone 
to detect and select supernova candidates, assigning them 
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galaxy photometric redshifts from an external catalogue. 
This method was applied in a deferred analysis of the first 
three years of the Supernova Legacy Survey (SNLS) con- 
ducted at the Canada-France-Hawaii Telescope (CFHT) 
from 2003 to 2008, as part of the Deep Synoptic Survey 
of the CFHT Legacy Survey (CFHT-LS). 

SNLS was defined with an optimised observing strat- 
egy. Four one-square degree fields were targeted through- 
out the 5 to 7 consecutive lunations where they remained 
visible. In every lunation, each visible field was repeatedly 
observed every 3 or 4 nights during dark time. This rolling- 
search mode allowed supernova light curves to be measured 
with very good time sampling in four broadband filters, 
den oted qm, '"m, iu and zm, similar to the SDSS filter 
set (jRegnault et al.1 . l2009t ) and spanning the wavelength 
range from 400 to 1000 nm. Th ese measurements w ere 
made with the MegaCam imager (jBoulade et al.L 120031 ). a 
1 square degree array of 36 CCDs with 340 million pixels 
in total. 

The standard SN la selection for the SNLS cosmological 
analy ses, based on real-time detection (Perrctt et al], 
|2010D and fo ll ow-up spectroscopy can be found 
in lAstier et al.l (20061) for the first year data sample 
and in Guv et al. ()2010[ ) for the 3-year sample. Spectral 



scopic identification, which has a lower limiting magnitude, 
are compared in Section [5.21 



da ta from the SNLS fo ll ow-up spectros c opy a r e described 



Howell et all dm 



5a) . iBronder et al.l . 

(|2008D . iBalland et all (I2009D and iWalker et al.l (I2010D . 

Photometric selections are commonly used in SNLS to 
define spectroscopic follow-up priori tisation from re a l-time 
(and thus partial) light curves (iSullivan et al.l . [2006h 
or to enlarge the real-time SN la sample in order t o 
measure SN la rates or properties (jSullivan et al.l . l2006bD . 
Photometric SN classification is also expected to become 
a key challenge for future supernova projects by lack of 
spectroscopic resources to confir m the large amount of 
supernovae that will be detected (|Kessler et al.l . |2010() . 

The photometric analysis described in this paper is in- 
dependent from the SNLS real-time detection pipeline and 
the ensuing cosmological analyses. It uses different data 
processing and detection, and applies a coherent set of cuts 
on 3-year light curves of all detected events, whether with or 
without spectroscopic information. Benefitting from more 
complete light curves and less stringent minimal brightness 
constraints than the real-time spectroscopic selection, the 
photometric analysis aims at going deeper in magnitude. 
As such, it can provide a cross-check to estimate possi- 
ble biases of the standard SNLS SN la selection for cos- 
mological analyses. It is also the basis of a photometric 
selection of non-type la supernovae that allowed the core- 
collapse su£ernovarateto be measured at a mean redshift 
z ~ 0.3 (jBazin et al.ll2009l ). Ultimately, it aims at exploring 
the feasibility of a cosmological analysis with a photometric 
SN la sample, which will be the subject of a future paper. 

The outline of the paper is as follows. Data process- 
ing, transient event detection and photometry are described 
in Section [2j The selection of SN-like transient events and 
their association with host galaxy photometric redshifts are 
presented in Section [3] The selection of SN la candidates 
is discussed in Section |4l The characteristics of the SN la 
candidates are presented in Section |S] where selected events 
with and without spectroscopic identification are compared 
in the same range of magnitude and found to have simi- 
lar properties. Finally, distance estimates from the whole 
photometric sample and from the subsample with spectro- 



(120081). lEllis et al.l 



2. Detection and photometry 

2.1. Transient event detection 

Data corresponding to SNLS observations from March 1*^*, 
2003 to September 2 1*'', 2006 were analysed, using images 
pre-processed at t he end of each lunation with the Elixir 
pipeline at CFHT (jMagnier fc Cuillandrj . 120041) which pro- 
vided flat-fielding and fringe subtraction. These data in- 
clude the presurvey period (up to June l'^*, 2003) which 
was used for the commissioning of MegaCam and thus cor- 
responds to unstable observing conditions. 

The astrometric solution for the images of each of the 
four SNLS fields (Dl D2 , D3 and D4) was co mputed 
with the TERAPIX Q calibration to ol SCAMP (iBertinl . 
I2006D and the USNO-Bl.O catalogue (|Monet et all . I2003D 
The im ages were then resa mpled with the TERAPIX tool 
SWarp (jBertin et al.l . l2002t ) according to the previously de- 
termined astrometric solution. A set of about 20 images 
with the best photometric quality (based on seeing and ab- 
sorption considerations) was selected for each field and each 
filter {gM, fM, im or zm) and co-added with SWarp in or- 
der to define reference images. The selected images were 
all chosen in the first or second season, depending on the 
field. The same reference images were used throughout the 
processing. 

Cosmic rays, defects due to dead pixels or induced 
by the resampling procedure were searched for in the re- 
sampled images and in the reference images. The content 
of the corresponding pixels was replaced by the estimate 
of the local sky background. An image subtraction pack- 
ag e, TRITON (iLeguillouL |2003[ ). based on the algorithm 
of lAlard fc Luptonl (jl998l) . was then applied on every re- 
sampled image. In this procedure, the convolution kernel 
and the sky background were determined independently on 
eight identical non-overlapping tiles paving each CCD of 
the mosaic. The number of tiles was chosen to optimise the 
precision of the subtractions, spatial variations of the kernel 
and sky background requiring small tiles, whereas the ker- 
nel determination requires a sufficient number of bright ob- 
jects and thus large tiles. To select bright objects to be used 
in the convolution kernel dete rmination for a giyen tile , the 
TERAPIX tool SExtractor (,Bertin et Arnoutd . I1996D was 
applied to both the reference and the current images, with 
a detection threshold at 2cr w.r.t. sky background. Bright 
objects present in both selections and neither saturated nor 
too close to a tile boundary were kept, which left about 100 
objects per tile for images at the average seeing of 0.7 arc- 
sec. The kernel and sky background were then fitted in 
59x59 pixel regions centred on each object. Adjacent ob- 
jects whose signal regions were overlapping by more than 
20% were discarded. Subtractions were considered as valid 
only if the integral of the convolution kernel was found to 
be above a filter-dependent threshold. The rate of lost sub- 
tractions (either failed or not valid) was around a few %. 

In each field, subtracted images in the im fiher were 
stacked for each lunation. On average, 27 images entered 
each stack. SExtractor was applied on each stack to con- 
struct catalogues of objects exhibiting a variable flux - 



see 
software 



http://www.astromatic.net/software] for the TERAPIX 
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whether positive or negative - in the «m filter in any luna- 
tion. At least 4 adjacent pixels with a signal of more than 
2.5(7 w.r.t. sky background were required to confirm an ob- 
ject. As the same set of reference images was kept through- 
out the whole processing, some supernovae may have part 
of their signal included in the references. In that case, the 
events might appear as missing flux on each of the images 
and could have been missed if searching only for a flux ex- 
cess. Considering negative measured fluxes when building 
the above catalogues allowed us to recover those events. 

All lunation catalogues were then merged to produce 
the final detection map. In this way, any object detected on 
several lunations gave only one detection, with a position 
averaged over all lunation stacks. 

When applied to the 3-year data sample of SNLS, the 
above detection pipeline produced a total of about 300,000 
detections, dominated mostly by saturated signals from 
bright objects which were not perfectly subtracted. The 
efficiency of the pipeline is discussed in the next section, 
while photometry is described in section 12.31 and further 
selections to discriminate between spurious detections and 
supernova signals are the subject of section |31 

2.2. Detection efficiency 

In order to check the performance of the above pipeline, 
we studied the detection efficiency with Monte-Carlo gen- 
erated artificial images produced for the Dl field in the im 
filter that we used for the detection. From this study we de- 
rived an efficiency model that allowed us to determine the 
efficiency for the other fields, D2, D3 and D4, which had 
slightly different mean observing conditions. This efficiency 
model was then coupled with a simulation of synthetic light 
curves (described in Section [23111) that we set up to study 
our photometric selection criteria. 

The Dl Monte-Carlo images were obtained by adding 
216,000 simulated supernovae to real images, thus naturally 
reproducing the observing conditions (sampling freq uen- 
cies a nd photometric quality) of the experiment (Ripoch^. 
1200 /t) . An external catalogue of host galaxies identified 
from deep image stacks of the CFHT-LS Deep Fields 
(lllbert et al.1 . 1200611 was used to randomly choose the posi- 
tion of each simulated supernova, following the observed 
surface brightness of the hosts, and to set the redshift, 
equal to the photometric redshift of the host. The red- 
shift was restricted to the range between 0.2 and 1.2. The 
supernovae were attributed a random stretch, colour and 
intrinsic dispersion term according to the distributions ob- 
served with spectroscopically identified SNe la. The time 
of maximum luminosity was drawn uniformly within each 
roughly six-month long season of observation. The i^vz-band 
light curve of each artificial s upernova was computed us- 
ing SALT2 ()Guv et al.l , I200i1) and the supernova flux at 
each observation date as deduced from the light curve was 
added to the corresponding images. The images were then 
processed with the same pipeline as the real images, setting 
the backgrounds and convolution kernels to those measured 
on the original images in order to avoid biases in the sub- 
traction procedure due to the presence of the numerous 
simulated supernovae. The detection efficiency was then de- 
fined as the fraction of simulated supernovae recovered in 
the detection map at the end of the processing. 

The detection efficiency for a supernova with date toi of 
maximum light in ij\f-band depends on its peak magnitude 



nioi, on the seeing and sky background during nearby ob- 
serving times tk, and on the relative epochs toi — tk- The 
Monte Carlo images allowed us to study the efficiency de- 
pendence on these variables, in order to build an efficiency 
model to extrapolate from Dl to the other fields. As an ex- 
ample, the average first season detection efficiency, e(moi) 
is illustrated in Figure [T] for Dl. The efficiency is rather 
uniform over the CCD mosaic, except for CCD 3 which 
failed to work during three out of the six dark time peri- 
ods of the first season. The efficiency is nearly magnitude- 
independent, emax ~ 0.97, out to moi = 23.5 in the SNLS 
magnitude system 0. This is followed with a steep de- 
cline at faint magnitudes reaching an efficiency of 0.50 at 
TOoi ~ 24.7. Very similar efficiencies were obtained in the 
other two seasons. Figure [T] also shows the result of our 
efficiency model that reproduces the dependence of the ef- 
ficiency as a function of the SN la peak magnitude. This 
model is described in Appendix A. 




Fig. 1. Detection efficiency e(moi) for the first season of 
Dl as a function of the generated peak magnitude in im- 
The model (blue curve) is compared with Monte-Carlo esti- 
mates (symbols), showing separately the efficiency for CCD 
3 (triangles), or averaged over the rest of the mosaic (dots). 




Fig. 2. Detection efficiency averaged over the entire CCD 
mosaic for fields Dl, D2, D3 and D4 as a function of the 
generated peak magnitude in im- 



^ The SNLS magnitude system uses star BD -1-17 4708 as a 
flux standard, which, as compared with the standard reference 
Vega, has the advantage of known Landolt magnitudes, precise 
SED measurements and colours close to the a vergae colours of 
the L andolt stars used in the SNLS calibration (iRegnault et all . 
120091 1 
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The efficiency model was applied to fields D2, D3 and 
D4, using the observation dates on these fields and the in- 
formation of seeing and background for each exposure. Each 
CCD was treated independently to reproduce exactly the 
possible failures, whether electronic or due to the process- 
ing of the real images. The detection efhciency computed 
for each field over the observing period between March 2003 
and September 2006 and averaged over the entire mosaic 
is shown in Figure [21 Fields D2 and D3 have a smaller effi- 
ciency than Dl and D4. This is caused by the lack of data 
due to bad weather in February 2004 and in March 2006, 
when only D2 and D3 were observable. 

2.3. Photometry and calibration 

For each of the 300,000 detections, four-filter light curves 
were built from individual subtracted images. TRITON dif- 
ferential photometry with Point Spread Function (PSF) fit- 
ting was applied, imposing the position found in the 
detection map (see Section [2?T|) . Fluxes were normalised by 
the integral of the convolution kernel determined in the sub- 
tractions, in order to express all fiuxes on a common scale, 
that of the reference image. For each exposure and CCD, 
the PSF was determined as follows. Bright events were first 
detected on each resampled image using SExtractor. The 
flux in the pixel of maximum content and the integrated 
flux of the event allowed saturated stars and galaxies to be 
identified and rejected. This left about 100 non saturated 
stars in each CCD, which were used to define one PSF per 
CCD by averaging the star profiles computed in a seeing- 
dependent box. 

The measured fluxes were calibrated using the S NLS 
tertiary standards published in lRegnault et al.l ()2009D . We 
measured the tertiary star fiuxes directly on current im- 
ages with the same photometry algorithm and PSF as 
for transient events and, as for the latter, the star fluxes 
were normalised by the convolution kernel integral deter- 
m ined in the subtraction s. Using the magnitudes published 
in lRegnault et all (|2009f) . the tertiary star fluxes were then 
converted into zero-points for each exposure and a weighted 
average of these zero-points was computed over each CCD, 
rejecting zero-points above 5ct w.r.t. the median value. The 
measured star fluxes were compared with fluxes derived 
from the star magnitudes by applying the average zero- 
points and the relative flux difference was fitted as a func- 
tion of flux (assuming a linear variation) . In the whole range 
of tertiary star fluxes, the fltted relative flux difference was 
found to remain below . 3% w hatever the field or band. 
Note that iRegnault et al.l (j2009t ) also provides fiat-fielding 
corrections due to residual radial non-uniformities of the 
imager photometric response. These corrections were ap- 
plied to all measured fluxes, prior to any further treatment 
(conversion into zero-point for stars or calibration for tran- 
sient events). 

The calibrated fiuxes were then filtered. Residual defects 
from bad pixels or cosmic rays, detections due to satellite 
trails or saturated stars were identified on subtracted im- 
ages as bright objects selected by SExtractor with a de- 
tection threshold at 3a w.r.t. sky background. Photometric 
points whose PSF overlapped the spatial extension of any of 
the above objects were discarded. Flux measurements ob- 
tained under bad seeing conditions (above 1.2") were also 
eliminated, as well as flux measurements in nights where all 



exposures were of bad quality 0. At this point, flux errors 
were renormalised so that after filtering the distribution of 
the pulls of the exposure fiuxes w.r.t. the night flux aver- 
age be of variance 1 in the test-sample of 278 light curves 
matched to spectroscopically identified SNe la in our detec- 
tions (see Section 15^ . Then, in any light curve, points with 
fiux errors either too low (w.r.t. the expected flux uncer- 
tainty due to sky background) or too high were discarded. 
Only exposures with a flux within 3(t from the night flux 
median value were retained and only nights with at least 
two such exposures were kept. On average over all recon- 
structed light curves, the fraction of measurements removed 
by the above cuts is at most 15%. 

Flux measurements corresponding to selected exposures 
within the same night were weighted by their errors and 
averaged to define mean light curves in each of the four 
filters. Finally, in order to account for a possible contribu- 
tion of variable objects in the reference images, a common 
baseline, defined as the mode of the flux distribution, was 
computed in each mean light curve and then withdrawn 
from the latter. 

Table 1. Performance of the photometry used in this work. 



Band 


Full sample 


Optimal conditions 




Median 


MAD 


Median 


MAD 


Qm 


0.006 


0.016 


-0.0001 


0.014 


tm 


0.018 


0.028 


0.006 


0.020 


hi 


0.015 


0.023 


0.001 


0.016 


zm 


0.009 


0.054 


-0.002 


0.047 



Notes. The table shows the difference mi — m2 of the SALT2 
fitted peak magnitudes obtained with the photometry of this 
paper (mi) and with the SNLS-3year photometry (m2) : median 
magnitude difference and median absolute dispersion (MAD) for 
all SNe la with good quality light curves (left) and for the subset 
detected in optimal conditions, as described in the text (right). 



The fiux estimation was checked with the set of 278 
spectroscopically identified SNe la previously mentioned. 
Peak magnitudes at the date of i?-band maximum light 
from our light curves were compared with th ose derived 
from the SNLS-3year light curves published in iGuv et all 
(j2010[ ) which correspond to a different image processing 
and an optimised photometry. In order to restrict to light 
curves of good quality, both sets were submitted to sam- 
pling and S/N cuts as described in iGuv et al.l (|2010( ). The 
selected light curves were then fitted with SALT2 in each 
filter fixing the date of -B-band maximum light, Xi and 
colour at the values obtained from a multi-band fit to the 
SNLS-3year hght curves (see Section [23)) . Only photomet- 
ric points present in both sets of light curves were consid- 
ered in the fits. 

Differences in the peak magnitudes fitted in each 
band are summarised in the left-hand column in Table [TJ 
Magnitudes from our light curves are fainter by at most 
0.02 mag with dispersions between 0.02 and 0.05 mag. Two 
effects explain these differences. Some of the SNe had SN 

^ grade C from the TERAPIX visual image qual- 
ity control as found in the TERAPIX Spica database, 
http: / /clix.iap.fr/steeringgroup /qf/ 1 
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light included in our reference images, which had no impact 
on their detection but affected their photometry. Second, 
even for SNe with no light in the reference images, the 
magnitude difference was observed to be larger when the 
SN was located near the top or bottom edges of the mo- 
saic. Restricting to SNe with no light in the reference im- 
ages and located far from the top/bottom mosaic borders, 
closer performance between the two photometry techniques 
are obtained, as shown in the right-hand column in Table [T] 
(column labelled 'optimal conditions'). 

The dependence of the magnitude difference with the 
SN location is likely to be due to our using the position 
found in the «m detection stack as the pivot of the PSF 
photometry. This is not optimal, especially on image edges 
due to distortions induced by alignment which are expected 
to be more important there. The average im position reso- 
lution provided by our detection pipeline was found to be 
0.41 pixel in the full sample of SNe la and 0.32 pixel for 
SNe in optimal cond itions, to be com pared with 0.28 pixel 
in the photometry of' Guv et al.l (|2010[ ) which performs a si- 
multaneous fit of SN la positions and fluxes^. Position mea- 
surement inac curacy lead s to un derestimated fluxes. Using 
appendix B of lGuv et all (|2010[ ). we found that resolutions 
of 0.41, 0.32 and 0.28 pixel would lead on average to Im 
biases of 0.011, 0.007 and 0.005 mag, respectively. These 
numbers reproduce the order of magnitude of the effect re- 
ported in Table [T] 

Summarising, when all SNe are considered, fluxes in this 
analysis are reconstructed with resolutions of a few % and 
an uncertainty below 2% on the absolute flux scale, which 
is accurate enough to set up a photometric selection. 

2.4. Catalogue of stars and galaxies 

The analysis described in this paper uses a catalogue of 
nearly 2 million stars and galaxies present in the four SNLS 
fi elds. It combines t wo sources. The first one, published 
in lllbert et al.l (|2006f ) . provides positions of stars and galax- 
ies, as well as photometric redshifts for the latter. The sec- 
ond was obtained by applying SExtractor on our reference 
images and contains measurements of the object size. The 
two catalogues were combined into one. Events common 
to both catalogues were classified as star or galaxy based 
on their classification in each of the catalogues. Ambiguous 
cases of star vs galaxy classification (e.g. at fainter magni- 
tudes) were solved with the help of flux and size measure- 
ments from our reference images. Altogether, the catalogue 
consists of 94% galaxies, 5% stars and only 1% ambiguous 
objects. 

The lllbert et"al] (|2006[ ) catalogue also provides photo- 
metric redshifts for more than 520,000 galaxies with an 
AB magnitude brighter than 25 in im- The photometric 
redshift method was applied to the CFHT-LS multi-colour 
data of the four deep fields of the survey. It was trained 
with around 3,000 VVDS spectra observed with VIMOS in 
one of the CFHT-LS deep fields, in both visible and near 
infrared bands. The method reaches a redshift resolution of 
0.037 with 3.7% of catastrophic redshift errors for a sample 



Note that this procedure is ideal when a limited number 
of objects is to be reconstructed but would be too time con- 
suming when dealing with all transient events detected by our 
photometric pipeline 



selected at im < 24. All galaxy photometric redshifts used 
in the following are from this catalogue. 

2.5. Light curve fitter 

Throughout this paper, we use the SALT2 package as a 
SN la ligh t curv e fitter. The version is that described in 
iGuv et all (l2010h which was trained on a larger data sam- 
ple and has a higher re solution co l our va riation law than 
the original version of IGuv et all (l200l . SALT2 models 
the mean evolution of the spectral energy distribution se- 
quence of SNe la. The model was trained on Branch-normal 
supernova data, both multi-band light curves and spec- 
tra, covering low and high redshifts up to 0.7. The result- 
ing flux model of a SN la at a given redshift is a func- 
tion of four parameters. Those are the i?-band global flux 
normalisation factor, Xq, the date of 5- band maximum 
light and two intrinsic rest-frame parameters, a stretch- 
related parameter, Xi and a colour parameter, defined as 
the {B — V) colour offset at the date of i?-band maximum 
light with respect to the average colour in the training sam- 

pie, C={B- V)r,^a.- <B-V>. 

In this paper, SALT2 was used either to produce syn- 
thetic SN la light curves assuming a given cosmology or to 
fit observed light curves under the assumption that they 
come from a SN la. In that case, the results are magni- 
tudes, either in the MegaCam filters or in the rest-frame 
_B-band. Except in Section [2731 these magnitudes are esti- 
mated from a global fit to photometric points in all filters 
that are expected to provide a significant contribution at 
the event redshift. 

3. Supernova selection 

First, we apply a set of selection criteria to the light curves 
of all candidate detections in order to reject spurious ob- 
jects and select for SN-like events. The SN selection criteria, 
described in the next section, were designed to be efficient 
for all types of supernovae. Spectroscopically identified su- 
pernovae present in our sample (see Section 13.21) and syn- 
thetic SN la light curves (see Section 13. 3p were used as 
qualitative guidelines to define the cuts. 

3.1. Selection criteria 

The selection proceeds in four steps. The light curves were 
first searched for a significant fiux variation in order to se- 
lect variable objects. We then checked that the main vari- 
ation in each curve had a shape consistent with that ex- 
pected from a SN event and that there was no other flux 
variation away from the SN-like variation. Events likely due 
to stars were then removed. Eventually, sampling require- 
ments were applied to ensure good quality light curves were 
selected. These four steps are detailed below. 

3.1.1. Search for flux variations 

In each light curve and filter, a search for significant flux 
variations was applied. Variations were due to start with 
a photometric point of positive flux and significance above 
la. They ended if a point of negative flux had a signifi- 
cance above Icr, if two successive points of positive fluxes 
had significances below Icr or at season ends, since super- 
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novae are expected to show only short-term variations last- 
ing typically over three consecutive lunations. To reduce the 
amount of spurious detections, the most significant varia- 
tions found in the im and tm light curves were required to 
contain each at least three points and to have their dates 
of maximum flux within 50 days from one another. This 
reduced the number of detections by about a factor 6. 

3.1.2. SN-like variation 

Most events at this stage were either spurious detections 
or objects varying on a longer term than supernovae. To 
be considered further, light curves were then required to 
have a shape compatible with that of SN-like events. To 
test this, the main variation in each filter was fitted with 
the phenomenological form : 



fit) 



(1) 



where k indicates the filter. Besides a constant, c , the fit 
parameters are , which sets the normalisation of the vari- 
able signal, Tj^ii and r^jge, which define the fall and rise 

times of the variation, respectively, and Iq which is related 



to the date of maximum by 



^^i'''fall/'''ri6 



1). Fits were independent in the four bands and run over 
the entire 3-year light curves. While form [1] has no partic- 
ular physical motivation, it is sufficiently general to fit the 
light curve shape of all types of supernovae. 

To reduce the contamination by long-term variable ob- 
jects or by random fluctuations, we first required that 
points away from the main variation were compatible with 
a constant flux. To do so, points outside the time interval of 
the main variation in Im were compared with the constant 
from the previous fit. Their significances w.r.t. that con- 
stant were added in quadrature over all four filters and the 
sum was normalised to the total number of points in the 
sum. This defined the following variable, hereafter referred 
as to off- variation ■ 



Xlff 



1 



EE 



where k indicates the filter and Fj represents the flux as 
measured at point j in the light curve with an error aj. A 
cut was applied on Xoff ^ function of the maximum flux 
observed in i^- As a result, about half of the detections 
resulted in light curves which, outside the main variation, 
were more irregular than those of SNe. 

The shape of the light curve was then tested for con- 
sistency with ([T]). This was done in the zm fllter only, as 
a significant signal is expected in this filter at all redshifts 
and for all types of SNe. A cut was applied on the difference 
in reduced between a fit by a constant and the fit by for- 
mula ([1}, as a function of the reduced of the latter fit. 
This cut is illustrated in Figure |3] where the events selected 
by all previous cuts are compared with simulated SNe la 
(this comparison is further discussed in Section [3.3.2p . Also 
highlighted in the plots are the spectroscopically identified 
events contained in our sample. A clear separation is ob- 
served between SN events and the bulk of the detections. 
After this cut, the number of events was reduced by a fac- 
tor ■~10. While this cut is more discriminating than that on 





Fig. 3. Difference in reduced between im light curve fits 
by a constant and by an SN-like shape (see text), as a func- 
tion of the latter x^- All cuts of the SN selection previous 
to that illustrated in this figure have been applied. Black 
dots stand for SNLS data (top) and synthetic SN la events 
(bottom). In both plots, green squares (resp. red triangles) 
are data events identified by spectroscopy as Type la (resp. 
core-collapse) SNe. Events above the curve are selected by 
the analysis. 



the off- variation x^, they proved to be complementary, as 
part 12%) of the events rejected by the first cut would 
pass the second one. 
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3.1.3. Star rejection 



Table 2. Effect of the SN selection cuts. 



Events at this point were still contaminated by stars, ei- 
ther truly variable stars or bright stars leading to saturated 
signals at some epochs. We discarded detections within 3 
ApsF from stars in the catalogue described in Section 
where ApsF is the half- width of the PSF in the im refer- 
ence images. Saturated signals from bright stars lead to 
light curves with large flux changes over a few epochs. 
Events whose zm flux varied by more than 70% between 
the date of maximum flux and the two closest epochs in 
the same season were thus rejected. The above two cuts re- 
duced the number of remaining detections by 12%. If they 
were applied instead of the constraint on the off- variation 
X^, the reduction factor at that level would be similar (46% 
instead of 49%) but the sample remaining after the selec- 
tion on the difference in reduced would increase by 65%. 
The three cuts previously described are thus necessary to 
achieve a significant rejection against background. 



3.1.4. Sampling requirements 

Quality criteria ensuring sufficient temporal sampling of the 
light curves complete the selection. In the Im and tm bands, 
at least one pre-max epoch and one post-max epoch were 
required in an interval ranging from 30 days before the 
date of maximum flux in each filter up to 60 days after 
that date. These conditions were found appropriate for SN 
events, even fast-declining ones like SNe la which on aver- 
age have time constants Trise ~ 4 days and Tfaii ~ 20 days 
in the im and tm bands. In the gM and zm bands, at least 
two epochs were required in a similar interval around the 
date of maximum flux in the im filter. Finally, the dates 
of maximum f^^, and i^^, from the fits by formula ^ 
were required to be within 50 days of each other, as more 
than 99% of both the synthetic SNe la and the spectroscop- 
ically identifed SN data events were observed to satisfy this 
condition. 



3.1.5. Results 

At the end of the SN selection, a set of 1483 events is re- 
tained. The effect of the cuts on data is detailed in Table [5J 
which also gives the output of the cuts on spectroscopically 
identified supernovae and on simulated SN la light curves 
(as described in the next two sections). 

3.2. Comparison with tiie SNLS real-time selection 

T he SNLS real- t ime S N selection procedure is described 
in iPerrett et al.l ()2010f ) . It relies on real-time detection of 
variable events based on photometry, combined with a spec- 
troscopic follow-up of the candidates which were both likely 
SNe la, as d etermined from an a nalysis of their real-time 
light curves (| Sullivan et al.l . l2006[l , and bright enough ones 
to allow for conclusive spectroscopic measurements. Events 
fainter than iff = 24.4 or far beyond maximum light at 
detection time were not observed spectroscopically. Events 



with moderate brightness, 22.9 < iff < 24.4 were quali- 
fied for spectroscopy if their fractional increase in brightness 
compared with their hos t galaxy was a bove a magnitude- 
dependent threshold (see iPerrett et al.l (|2010f ) for more de- 
tails). The spectroscopic observations were essential to con- 
firm the nature of the real-time candidates and to obtain 



Cut 


Events 


la 


CC 


e(Ia) 


Detection in im 


295683 


278 


55 


0.96 


Variations in zm and tm 


193104 


277 


55 


0.94 


Consistent tmax in im and 


52284 


276 


53 


0.93 


Cut in Xo// 


26510 


276 


52 


0.93 


2 2 ■ ■ 

Xconstant " Xshapc Cut in IM 


2474 


276 


52 


0.93 


Star veto 


2165 


275 


52 


0.92 


Sampling cuts in jm, tm 


1674 


253 


43 


0.79 


Sampling cuts in gM, zm 


1571 


246 


42 


0.73 


Consistent t^'^a; in ^a/ and tm 


1483 


246 


42 


0.73 



Notes. The table details the effect of the cuts on the first three 
years of SNLS data : all detections (column 2), subsamples of 
events identified by spectroscopy as Type la and core-collapse 
SNe (columns 3 and 4) . The last column indicates the efficiency 
of the cuts for bright SNe la, derived from synthetic SN la light 
curves at low magnitudes (moi < 23). 



precise redshifts. As an example of how efficient the selec- 
tion of targets for spectroscopy was, 71% of the first-year 
high redshift events sent for spectroscopy at the Gemini 
telescopes were confirmed to be Type la supernovae, either 
secure (" SN la") or probable la events ("SN la?"), as de- 
scribed in lHowell etall ([2001). 

In this section, we compare our photometrically iden- 
tified SN-like events to those located by SNLS in real- 
time. The SNLS databas^B lists all variable events tar- 
geted by the real-time pipeline, as well as information re- 
turned by spectroscopy. Snapshots of these databases were 
extracted in January, 2009 and used to match the se- 
lected SN-like events of our photometric analysis with SNLS 
real-time candidates. Spectral data u sed for this matching 
(Howeh et al.l (l2005l). iBronder et al.l (|2008[ ) . iBalland et al.l 
(|2009D andlWalkcr et al.' ('2010^') are identical to that used 
to d efine the SNLS 3-ycar sample for co smological analy- 
ses (|Guv et all . I20T0I: iConlev etall . [20T1I) . 

Of the 1483 selected SN-like events, 865 (58%) could 
be associated with real-time candidates. Most of the re- 
maining 618 events are faint events, with im peak magni- 
tudes above 23.4, as shown in FigurejH Besides background 
events, which are still present at this stage, these may also 
be interesting events too faint to be detected by the real- 
time pipeline. The matched events have the following types 
in the SNLS databases. 246 events are classified as SNe la 
and 42 are classified as core-collapse SNe, adding in both 
cases secure and probable events. These two subsamples are 
referred to as spectroscopically identified events throughout 
this paper. The other events are mainly classified as SNe 
(492 potential ones and 52 confirmed as supernovae by spec- 
troscopy). A few AGN's (12 potential ones and 3 confirmed 
events) and possible stars (11 events, none confirmed) are 
also present. The type of these events is mostly based on 
real-time photometry, as most of them were not sent for 
spectroscopy. The break down of the SNLS database event 
types in the selected SN sample is presented in Table [31 

Spectroscopically identified supernovae can be used to 
check the efficiency of the SN selection on bright supernova 
events. The SNLS database lists a total of 280 spectroscop- 
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Fitted peak i,^ magnitude 



Fig. 4. Fitted iM peak magnitude distribution of selected 
SN-like events. Among these, events matched to SNLS real- 
time candidates (resp. spectroscopically identified SNe la) 
correspond to the red (resp. green) hatched histogram. 



ically identified SNe la and 55 spectroscopically identified 
core-collapse events in the data covered by our photometric 
analysis. All of them are detected by our pipeline, except 
for two SNe la in D3, one outside our reference image and 
one from the presurvey phase. Of the 333 detected events, 
288 passed the SN selection, as shown in Table [5] Most of 
the lost events were rejected by the sampling cuts. Those 
were not fulfilled due to either bad weather conditions or 
observations during the presurvey phase (42%), or trunca- 
tion of the light curve at one end of the observation season 
(29%) or photometric points missing near maximum light 
because of CCD failures or full moon (24%). 

3.3. SN selection efficiency for Type la supernovae 

The efficiency of the SN selection was studied further with 
our set of synthetic SN la light curves. Those were simu- 
lated in the observed bands of MegaCam by the SALT2 
package. Details about the simulation are given in the sec- 
tion below, comparison with data and selection efficiency 
are discussed in the next two sections. 

3.3.1. SN la light curve simulation 

Synthetic SN la light curves were produced assuming a flat 
ACDM cosmology with Hm = 0.23. Redshifts were simu- 
lated in the range from to 1.2, with distances computed 
in the above cosmology. A volumetric distribution was as- 
sumed in order to reproduce the acceptance of the survey. 
The rcdshift dependence of the SN explosion rate is irrel- 
evant for efficiency studies and was thus ignored at this 
stage. Values of the SALT2 Xi and colour parameters were 
generated according to model distributions fitted to the 246 
spectroscopically identified SNe la present in our selection 
(see Table [5]), restricting to redshifts below 0.7 to avoid 
selection biases. A Gaussian distribution was used for the 
colour parameter, with mean and a values of 0. and 0.1 
(see Figure O. For Xi, a double Gaussian was found more 



appropriate, with mean values of 0.06 and -1.07, a values 
of 0.8 and 0.7 and relative amplitudes equal to 0.25. 




Fig. 5. Xi and colour distributions of spectroscopically 
identified SNe la with rcdshift z < 0.7. The curves are 
used to produce synthetic light curves from the observed 
distributions. 



For each generated event, a rest- frame i?-band peak 
magnitude, m*g, was computed as a function of the gener- 
ated values of Xi, C and of the luminosity distance at the 
generated redshift in the assumed cosmology, dL{z, J^m) as: 

m*B = M - aXi +(3C + 5\ogwdLiz, VIm) 

using M = -19.09, a = 0.13 and /3 = 2.56. A Gaussian 
spread with a — 0.14 was applied to m*g in order to al- 
low for SN la intrinsic dispersion. The magnitude after 
dispersion was converted into a global normalisation fac- 
tor Xq according to the relationship between m*g, Xq, C 
and Xi measured in the spectroscopically identified SNe la. 
Finally, dates of i?-band maximum light were drawn uni- 
formly in intervals starting 10 days before the beginning of 
the true observation periods and ending 5 days after their 
end. Each generated SN was also attributed a random po- 
sition in one of the four fields of observations. Corrections 
for the Milky Way extinction of the SN la flux at the 
generated co ordinates were applie d in SALT2 using the 
dust maps of ISchlee el et al.l (jl998|) and the extinction law 
of[Cardelli et^ (1989(1 

Detection and instrumental effects were then simulated 
in the following way. The detection efficiency model de- 
scribed in Section [52] was used to compute the detection 
probability of each synthetic event. Based on this proba- 
bility, generated events were then randomly eliminated or 
kept for further processing. For each event passing this step, 
fluxes were computed at dates corresponding to true obser- 
vations at the generated SN position in the mosaic. This 
allowed us to account for the observing strategy as well 
as for periods of lack of observations due to bad weather 
conditions or temporary CCD inefficiencies. 

From these ideal fluxes, reconstructed fluxes and errors 
were then simulated according to reference distributions de- 
termined from the 278 spectroscopically identified SNe la 
that we detected (see Table [5]). Distributions were those of 
flux errors and differences between the measured and the 
SALT2-estimated fluxes (thereafter called shifts in flux). 
Correlations between flux errors, flux values and flux shifts 
were taken into account when building the reference dis- 
tributions. As around 20% of the light curves in the test- 
sample exhibit more points with a large shift in flux than 
expected from a pure random distribution, the reference 
shift distributions were determined separately for good and 
bad quality light curves. Observation dates showing system- 
atic large shifts and errors were also described by reference 
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distributions determined separately. A total of 20,000 syn- 
thetic light curves was generated. 

3.3.2. Comparison with data 

Synthetic light curves are compared with spectroscopically 
identified SN la data in Figure [51 where the first step of the 
selections described in Section lOl and an observed peak im- 
band magnitude below 23.5 were required. The agreement 
is reasonable. In particular, the simulation reproduces both 
the peak and the tail of the distributions of the various 
variables used in the selections. 
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Fig. 6. Comparison between spectroscopically identified 
SN la data (open histogram) and synthetic light curves 
of bright SNe la (hatched histogram) in the Dl field. 
Distributions show the total number of points in the iu 
light curve and three variables used in the analysis (see 
text). The mean and RMS of each distributions are indi- 
cated in the upper box for for data and in the lower one for 
simulation. 



The comparison is extended to all data in Figure |31 
showing the qualitative agreement between synthetic SN la 
light curves and identified SN data. In particular, there are 
two features that are reproduced by the simulation. The 

of the analytical fit, X/it, deteriorates at large x^ differ- 
ences, that is at low redshift, where formula ([T]) is less suited 
to fit SN light curves due to a secondary maximum in the 
near-infrared rest-frame light curves that SNLS observes in 
the ij\/-band at low redshift. Second, at moderate Ax^, the 
X^jj distribution has a tail at values above 2 which, in the 
simulation, is reproduced by bad quality light curves, i.e. 
light curves with a sizeable number of photometric points 
departing from the generated fiuxes (see Section fS.S.ip . 

3.3.3. Efficiency 

The selection efficiency for Type la supernovae is presented 
in Figure [7] as a function of the generated peak magnitude 
in iM-band (full dots). The four fields were combined and, 
as fiducial cuts, we restricted all efficiency computations 
to those events whose generated date of _B-band maximum 
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Fig. 7. Efficiency of the SN selection for Type la super- 
novae as a function of the peak magnitude in jj\f-band, for 
all events (full dots) and for events with well sampled light 
curves (open dots). The efficiency was evaluated with syn- 
thetic SN la light curves. 



light was strictly within the observing periods. For com- 
parison. Figure [7] gives also the efficiency for well sampled 
light curves (i.e. events passing the sampling cuts of the 
selection), showing that, at this stage of the analysis, the 
efficiency loss for bright events is due mostly to sampling 
requirements. 

The efficiency for bright synthetic events is detailed in 
the last column of Table which also gives that for spec- 
troscopically identified SNe la in data. The efficiency from 
simulation (0.73/0.96 = 76%) is lower than that for data 
events (246/278 = 88%) since these have been submitted 
to tighter quality criteria to ensure conclusive spectroscopic 
measurements. As an example, bright events with a light 
curve starting past maximum light were not sent for spec- 
troscopy, while they are included in the simulated sam- 
ple. Restricting again the comparison to well sampled light 
curves in both simulation and data, the efficiencies become 
246/246 = 100% in data and 0.97/0.98 = 99% in simula- 
tion, which gives confidence that reliable selection efficien- 
cies are deduced from the simulation. 



3.4. Redshift assignment 

The selected SN-like events were matched with galaxies in 
the catalogue described in Section The match was con- 
sidered successful if the nearest galaxy from an event was 
found within a distance of Srgai, where the galaxy's effective 
radius, rgai, was defined as the half- width of the galaxy in 
the direction of the event and derived from the A,B and 9 
SExtractor parameters. In this procedure, galaxies so faint 
that no size measurement could be performed (typically 
fainter than (m ~ 24.5) were assigned a size of 1 pixel 
and were considered only if no bright galaxy was found in 
the vicinity of the event. The choice of the 57'gai threshold 
was a compromise between host finding efficiency and ac- 
cidental mismatching. Of the 1483 selected SN-like events, 
1352 (91%) had matched hosts and 1233 of these (91%) had 
a photometric redshift which was then attributed to the 
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event. The types of these events as found in the real-time 
database of SNLS (see Section 13. 2p are listed in Table [3l 

Host galaxy photometric redshifts were compared to 
event spectroscopic redshifts in the subsample of 346 se- 
lected SN-like events which have both measurements. The 
result is shown in Figure [H] In most cases, the two red- 
shift measurements agree and the galaxy photometric red- 
shifts provide a resolution of o'az/(i+z) ^ 3% in the cen- 
tral part of the distribution. However, 5.2% (resp. 10%) of 
the galaxy photometric redshifts depart from the spectro- 
scopic redshifts by more than 5 (resp. 3) a. Such differences 
can be explained by incorrect SN-to-galaxy associations or 
by bad galaxy photometric redshifts. The latter effect is 
li kely to be th e main reason for the deviations. According 
to lllbert et al.l (j2006[ ) 3.7% of the galaxy photometric mea- 
surements are to be expected with Az/ (\ + z) > 0.15, for a 
sample of galaxies with magnitudes below 24 in iff , which 
corresponds to the range of host galaxy magnitudes in this 
study. 




-0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 



Az/(1+z) 

Fig. 8. Host galaxy photometric redshifts compared with 
event spectroscopic redshifts in the subsample of events 
with both measurements available at the end of the SN se- 
lection. The resolution of galaxy photometric redshift mea- 
surements is given in the bottom plot. 



Incorrect redshift assignments are likely to generate sig- 
nificant systematic biases in the subsequent photometric 
SN la selection, especially at large {z > 0.6) redshifts, where 
most of the incorrect assignments occur (see Figure [8]). 
Special care was thus devoted to identify and reject most 
of these events in the SN la selection criteria that are de- 
scribed in the next section. As synthetic light curves provide 
an important tool to define the selections, host galaxy pho- 
tometric redshifts were also simulated for the sample of syn- 
thetic light curves. Redshifts were generated according to a 
Gaussian distribution in Az/{l + z)'^, to account for the fact 
that the reliability of the galaxy photometric redshifts de- 



creases significantly for fainter objects (lllbert et al.l . l2006h . 
Redshift outliers were generated with a Gaussian distribu- 
tion in Az/(1 -|- z) and assuming their rate rises linearly 
with redshift. Parameters for this simulation were derived 
from the comparison between host galaxy photometric red- 
shifts and event spectroscopic redshifts in the data, using 
the test-sample of 346 events already mentioned. Finally, 
the fraction of simulated events with a host redshift after 
the SN selection was given the same probability as observed 
in the data. 

Table 3. Event type break down at different analysis 
stages. 



Event 


Analysis step 




cIeiss 


SN sel. 


+Zgai la sel. 


All events 


1483 


1233 


485 


In database 


865 


739 


388 


la type 


246 


208 


175 


CC types 


42 


38 


1* 


SN-like 


544 


467 


210 


AGN-like 


15 


13 


2 


Variable-like 


11 


8 





Unknown type 


7 


5 






Notes. The break down is given in the photometrically selected 
SN-like events (column 2), in those events with a host galaxy 
photometric redshift (column 3) and in those finally kept as 
SN la events (column 4). The type of the event marked with an 
asterisk was revisited after this analysis (see Section |4]2]). 



3.5. Summary 

The total sample of 1233 SN-like events with a host galaxy 
photometric redshift is the starting point of the SN la 
search described in the next section. A measurement of 
the rate of core-collapse supernovae at low redshift was 
also d educed from that sample, as described in lBazin et al.l 

(|2nn9h .FI 

4. SN la selection 

The photometric SN la selection cuts were designed to dis- 
criminate SNe la from core-collapse supernovae (SNe CC) 
and to assess the reliability of the redshift assignment. The 
search was restricted to normal SNe la only, and relied on 
the SALT2 light curve fitter. As for the SN selection, the 
criteria were set up using spectroscopically identified su- 
pernovae present in our sample and synthetic SN la light 
curves as qualitative guidelines. 

4.1. Selection criteria 

The four-band light curves of each event were fit simulta- 
neously with SALT2, fixing the redshift to the host galaxy 
photometric redshift, to derive the date of B-band maxi- 
mum light i^o^, the colour C, the stretch-related param- 
eter Xi and the rest-frame i?-band peak magnitude m|j 

" This paper used a preliminary flux calibration resulting in 
slight difi^erences in the analysis (e.g. 1207 SN-like events with 
a host galaxy photometric redshift instead of 1233) 
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under the assumption that the event was a SN la. The fit 
converged for 98% of the events. 

4.1.1. Light curve sampling 

To ensure that meaningful SALT2 parameters were ob- 
tained, tighter quality cuts were applied w.r.t. the SN se- 
lection. Defining the event rest-frame phase as r = (io&s — 
t^ax)/{^ + ^gai), events were required to fulfil the following 
minimum conditions Q : 

(i) at least one measurement in the range — 10 < r < -1-5 
days for a reasonable estimate of tmax and thus of the 
peak magnitude, 

(ii) at least one measurement in the range -1-5 < t < -1-20 
days for a reasonable shape evaluation, 

(iii) at least one colour among {g — i), (r — z) and {i — z) 
with at least one measurement in each band in the range 
— 10 < T < -1-35 for a reasonable la/CC discrimination 
(see Section liXi)) . 

After these cuts, the selected sample contained 1152 events, 
among which 203 (resp. 35) were confirmed as SNe la (resp. 
SNe CC) by spectroscopy. The remaining selections in- 
tended to enrich the sample in actual SN la events with 
correct assigned host galaxy photometric redshifts. 

4.1.2. requirements 

To reject SN la light curves with a poor fit and to start 
eliminating non-la SNe, cuts were applied to the per de- 
gree of freedom in each filter, < 10 in the g filter and 

< 8 in the r, i, z filters. In addition, the x^ P^r degree 
of freedom of the overall fit was required to be less than 6. 
Seven spectroscopically confirmed SNe la were rejected by 
these conditions. Two were assigned an incorrect host pho- 
tometric redshift, the other five exhibited noisy light curves. 
These cuts also discarded nine spectroscopically confirmed 
SNe CC, all of the SN II type. 

The light curve fit with SALT2 only included MegaCam 
filters corresponding to the redshiftcd wavelength interval 
between 290 nm and 700 nm, wher e the SN l a mod el is 
defined with sufficient accuracy (see iGuv et al.l (l2007l^ for 
details). Filter gu was thus only used below a redshift of 
about 0.68, while zm was considered only for a redshift 
larger than 0.26. Outside these redshift intervals, SNe la are 
expected to give no significant signal in these bands, but 
other types of supernovae may have a different behaviour. 
Comparing the observed fiux in the unfitted bands to that 
corresponding to the SALT2 fit can give a complementary 
constraint against core-collapse events, as well as against 
incorrect redshift assignments in true SNe la. 

The x^ per degree of freedom with respect to the SALT2 
best fit was thus computed in the unfitted bands and re- 
quired to be less than 6 in g]\[, and less than 3 in zm- This 
criterion is stronger than that for the bands entering the fit, 
to compensate for the fact that the subsequent selections 
are looser for events which do not have all bands included in 
the fits. Ten spectroscopically confirmed SN la were elimi- 
nated by this second set of cuts. They were all assigned a 
bad host photometric redshift, with {\zgai — Zspe\) — 0.35 
compared to an average of 0.06 for the other confirmed 



SNe la. These cuts also rejected four spectroscopically con- 
firmed core-collapse events (two SNe II, two SNe Ib/c) be- 
cause of their high x^ in zm- 

4.1.3. Stretch and colour constraints 

As a further discrimination between SN la and non SN la 
events, a constraint was applied on the fitted Xi and colour 
parameters, requiring 



-0.2 

lo 



C 
(135 



< 1 



as illustrated in Figure [S] The constraint accounts for the 
offset of 0.2 observed in the stretch distribution, as shown 
in Figure [5] Such an offset may reflect the higher num- 
ber of high stretch SNe la observed at high redshift due 
to the increased fraction of star forming host galaxies in 
the pa st universe, which lead to more luminous super- 
novae (|Howell et all . 120061: ISullivan et all . I2010D . 

While the lower bound on Xi only removed a few rela- 
tively short and faint (average peak magnitude in im of 25) 
data events, the upper bound discarded a large number of 
long duration events among which seven spectroscopically 
confirmed SNe CC, all of SN II type. On the other hand, 
the constraint on Xi removed a negligible fraction (1%) of 
the synthetic SNe la. 
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These requirements are very simila r to those used to define 
the SNLS 3-year cosmological samples (|Guv et al.l . 120101 ) 



Fig. 9. Xi vs. C for synthetic SNe la (blue dots) and for 
data (all other symbols), after sampling and x^ constraints 
have been applied. Green filled circles and red triangles 
stand for spectroscopically identified SNe la and SNe CC 
present in our sample, respectively. Open black circles stand 
for data events with no spectroscopic identification. Events 
outside the red ellipse are rejected. The pileup at Xi=5 is 
due to an upper bound encoded in the fitter. 



The constraint on the fitted colour relies primarily on 
the observation that, while SN la colours were centred 
around (with a r.m.s. of 0.11), the spectroscopically iden- 
tified core-collapse events in our sample exhibited an aver- 
age colour (C) = 0.3. Requiring C < 0.35 rejected nine of 
these, mostly SN lb and SN Ic events. 

The cut on colour also allowed the rejection of SNe la 
where the associated host galaxy photometric redshift, Zgai, 
was far from the true one. As an illustration. Fig. [TU] 
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presents the difference between Zgai and the more pre- 
cise spectroscopic redshift, z^pe as a function of the fit- 
ted event colour, for events with both redshift assignments. 
Extreme colours in identified SNe la are mostly associated 
with inaccurate galaxy photometric redshift measurements. 
Synthetic SNe la show the same trend. The accuracy with 
which the colour of a true SN la is recovered, is indeed di- 
rectly related to the reliability of the redshift assignment. 
As an example, about 70% of the synthetic light curves with 
\Zgal ~ Zspe I > 0.2 have |C| above 0.35. 
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Fig. 10. Difference between redshift assignments as a func- 
tion of SALT2 fitted colour. Same colour code and level of 
selections as in Fig. [HI except for pink crosses which indi- 
cate data events rejected by the selections in unfitted 
bands. Red lines indicate the extreme values of the colour 
cuts. 



4.1.4. SALT2 colour-magnitude diagrams 

The last criteria to reject non SN la events were based on 
colour-magnitude diagrams from SALT2 fitted magnitudes. 
In these diagrams, Type la supernovae populate a thin band 
while core-collapse supernovae lie in a broad region which 
is shifted w.r.t. the SN la band. To emphasize the colour 
and magnitude difference between the two classes of events, 
the following diagrams were chosen: g — i vs. g, r — z vs. r 
and i — z vs. z (Fig. 1111 top to bottom), where g, r, i and z 
are the magnitudes in the MegaCam filters computed from 
the SALT2 best fit model, at the date of B-band maximum 
light. Note that events entered these diagrams only when 
the filters of interest were actually considered in the fit (see 
Sec. l4.1.2| . The first diagram was thus used for low redshift 
events, 0.16 < z < 0.68, the second one was considered in 
the redshift range 0.26 < z < 1.15, while the third one 
applied for 0.26 < z. Events at the highest redshift {z > 
1.15) were thus tested only through the third diagram. 

The g — i vs. g diagram proved to be very helpful in 
reducing the remaining background from core-collapse su- 
pernovae, as most of those we could detect were at low 
redshift, and, in this diagram, appeared fainter than actual 
SNe la for the same g — i colour. Although the magnitude 
dispersion of the SNe CC is close to 1 mag, this diagram 
offered a good discrimination between CC and Type la su- 




fitted i 



Fig. 11. Colour- magnitude diagrams based on SALT2 fit- 
ted magnitudes for events that passed all previous con- 
straints up to Xi and C cuts. Same colour code as in Fig.[9l 
Events below the lines were rejected, la magnitude ranges 
are indicated for 0.1 redshift bins centred on the indicated 
redshift values (vertical locations of the segments are arbi- 
trary). 



pernovae. Part of the SNe CC were detected at high enough 
redshift to enter the r — z vs. r or i — z vs. i diagrams. A 
similar trend was observed there, with a population dis- 
tinct from the SNe la, lying at fainter magnitudes for a 
given colour. The separation between the two populations 
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in these diagrams is however less pronounced than in the 
g ~ i vs. g diagram. 

Altogether, the three diagrams were useful to complete 
the reduction of the CC supernova background. Note how- 
ever that the spectroscopically confirmed SNe CC that 
would be rejected by the colour-magnitude diagrams alone, 
i.e. not considering any of the previous SN la selection cri- 
teria, are all typed as either SN-II or SN-IIP. In the colour- 
magnitude diagrams, the SNe Ib/c often lay close to the 
upper limit of the SN la band and could not be identified 
unambiguously. 

4.1.5. Results 

At the end of the selection, 485 SN la candidates were se- 
lected. The overall selection efficiency is given in Figure \T2\ 
as a function of the peak magnitude in the iM filter, for all 
synthetic events (full squares) and for events with well sam- 
pled light curves (open squares) defined as in Section [5.3.3l 
For the latter, the bright event efficiency is 67%. It would be 
80% if all events were assigned a host galaxy photometric 
redshift. This number represents the intrinsic performance 
of our selection for bright SNe la, regardless of the survey 
observing strategy or the availability of host redshifts. 
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Fig. 12. Efficiency of the SN la selection for Type la su- 
pernovae as a function of the peak magnitude in iM-hand, 
for all events (full squares) and for events with well sam- 
pled light curves, when the host galaxy redshift assignment 
probability is included (open squares) or not (crosses). The 
efficiency was evaluated with synthetic SN la light curves. 



Table 4. Effect of the SN la selection cuts. 



Put 


JZj VcIILd 


lei 




e^iaj 


SN selection 


1483 


246 


42 


0.97 


Host redshift available 


1233 


208 


38 


0.80 


Light curve sampling 


1152 


203 


35 


0.75 


x' 


951 


186 


22 


0.68 


Xi and colour 


596 


176 


7 


0.67 


g-ivs.g 


539 


176 


1* 


0.67 


r — z vs. r 


490 


175 


1* 


0.67 


i — z vs. z 


485 


175 


1* 


0.67 



Notes. The effect of the cuts is given for all detections (column 
2), for the subsamples of events identified by spectroscopy as 
Type la and core-collapse SNe (columns 3 and 4), and for syn- 
thetic bright SNe la (moi < 23) with well sampled light curves 
(column 5). The type of the event marked with an asterisk was 
revisited after this analysis (see Section [4. 2|l . 



To check this, we show in Figure [13] the selection efficiency 
as a function of the generated redshift (see Appendix B 
for the colour and Xi variations). The dependence seen in 
Figure [T3l follows mostly from the fact that the selection is 
easier for brighter SNe la, i.e. for lower-redshift, bluer and 
higher-stretch events. The SN la selection efficiency starts 
to decrease at redshifts above 0.6 while the SN selection 
efficiency plateau extended up to z ~ 0.8. We checked that 
this decrease is due to the cuts on colour and on the re- 
duced from the unfitted bands, which are essential to 
reject non-SNIa contaminants and true SNe la which were 
assigned an incorrect redshift. 
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The effect of the different selection criteria is detailed 
in Table |H where the selection efficiency is given for bright 
events with well sampled light curves. As for the SN selec- 
tion, the efficiency of the SN la selection criteria from sim- 
ulation (0.67/0.97 — 69%) agrees well with that for spec- 
troscopically identified SNe la (175/246 = 71 ± 3%). 

4.1.6. Discussion 

The use of colour-magnitude diagrams in the selection may 
generate biases in the selection as a function of redshift. 



Fig. 13. Selection efficiency from synthetic SN la light 
curves as a function of the generated redshift, at different 
stages of the analysis. 



The colour-magnitude diagrams used in the above anal- 
ysis are based on SALT2 fitted magnitudes and not on 
colours and magnitudes derived directly from observations. 
We checked that using magnitudes at the SALT2 date of B- 
band maximum light derived from fits to the light curves 
with formula (H]) would indeed lead to a poorer discrim- 
ination between SNe la and contaminants, due to larger 
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dispersions, especially in the r — z vs. r and i — z vs. i 
diagrams. More quantitatively, applying the selections de- 
scribed in Section 14.1.41 on these diagrams would lead to 
6% less efficiency and 75% more contamination by CC SNe 
(see Section 133]). 
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Fig. 14. Hubble diagram for data events at difi^erent levels 
of the photometric SN la selection. Black dots stand for 
events at the end of the selection, green filled (resp. open 
blue) circles for events rejected by the constraints on the 
SALT2 fitted g — i vs g (resp. i — z ys i) colour-magnitude 
diagram, open red circles for events rejected by the con- 
straints on the SALT2 fitted r — z vs r diagram. 
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Fig. 15. Host galaxy photometric redshifts compared with 
event spectroscopic redshifts in the subsamplc of events 
with both measurements available after all SN la selec- 
tion criteria have been applied (black dots) . Downward light 
blue (resp. upward pink) triangles stand for SN-like events 
rejected by the (resp. colour) cuts of the SN la selection. 
The redshift resolution at the end of the SN la selection is 
given in the bottom plot. 



The impact of the colour-magnitude constraints based 
on SALT2 fitted magnitudes is illustrated in Figure [T3] 
which shows the Hubble diagram of the sample of pho- 
tometrically selected events (see Section [5] for more de- 
tails about the distance modulus computation) as well as 
events rejected by each constraint. On average, the aver- 
age distance modulus of the rejected events is 1.34 mag 
larger than that of the selected sample, with a r.m.s. of 
0.61 mag, almost twice that of the selected sample. Keeping 
events whose distance modulus is within 0.55 mag above the 
Hubble fiow would lead to similar SN la efficiency, core- 
collapse contamination and sample of data events as the 
colour-magnitude diagram constraints used in this analy- 
sis. The two methods are thus equivalent. 

As mentioned in Section 13.41 redshift outliers were 
present in a significant fraction of the selected SN events, 
e.g. 10% of those had redshifts which were 3 a outliers. 
Figure fTSl shows the comparison between spectroscopic and 
host galaxy photometric redshifts at the end of the SN la 
selection. Highlighted in the plot is the effect of the and 
colour requirements, which prove to be helpful to reject 
events which were assigned a redshift outlier. At the end of 
the SN la selection, 1.4% (resp. 5.5%) of the host galaxy 
photometric redshifts depart from the spectroscopic ones 
by more than 5 (resp. 3) cr. The outlier fractions at the SN 
selection level have been reduced by a factor 4 and 2, re- 
spectively. There is no change in the redshift resolution, but 
the redshift distribution is now much closer to a Gaussian, 
as shown in the bottom plot of Figure [TSl 



Summarising, with the above photometric selection, 
well sampled light curves of bright SNe la are selected 
with an average efficiency of 67%, which includes the 83% 
efficiency of the host galaxy photometric redshift assign- 
ment. Photometric redshifts derived from supernova light 
curves would offer a better re dshift assignment probability 
than w hat we achieved here. In lPalanaue-Delabrouille et al.l 
we set up such a method based on SALT2 and tested 
it with the synthetic SN la light curves and the photometric 
data sample described in the previous sections. Combining 
both supernova selection and redshift determination from 
photometry will be the subject of another paper. 

In the next sections, we pursue the discussion on the 
present analysis, comparing its output to that of the real- 
time SNLS analysis and evaluating the contamination of 
the selected sample by residual core-collapse supernovae. 

4.2. Comparison with the SNLS real-time selection 

As in Section [3?2l events of the photometric sample were as- 
sociated with SNLS real-time candidates. The break-down 
of the event types is summarised in Table [3] Out of the 485 
candidates, 388 (80%) were detected in real-time. Of these, 
45% are spectroscopically identified SNe la (adding up se- 
cure and probable SNe la) and another 54% were typed 
(mostly photometrically) as SNe. The remaining 3 events 
have the following characteristics. 

One event, SNLS 06D2bo was identified as having a 
spectroscopic redshift of z=0.370, on the basis of strong 
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[Oil] and [OIII] emission lines from the host galaxy. The 
best match appeared to be the broad-lined SN 1998bw, so 
the event was classified as a possible CC SN. The event 
is highlighted in Figure [TT] (middle plot) as the red trian- 
gle at tm ^ 25. It lands within the locus of points de- 
fined by SNe la. The photometric redshift of the host was 
much higher, z = 0.82, and agreed with the estimate of 
the redshift derive d from fitti ng the SN light curve (using 
e.g. Sullivan et al] ()2006l ) or iPalangue-Delabrouille et al] 
(|201ClD ). The large discrepancy between the photometric 
redshifts (SN and host) and the spectroscopic one prompted 
us to re-examine the spectroscopic data. We soon discov- 
ered that the spectrum of an unrelated object had been 
extracted instead of SNLS 06D2bo. The data were repro- 
cessed. SNLS 06D2bo is clearly a SN la, with a best fit 
redshift of 2 = 0.79. We have chosen to let the error in the 
extraction stand. It vindicates the accuracy of the tech- 
nique we have developed to photometrically identify SNe 
la, and it indicates the level of human error in the analysis 
of the SNLS spectroscopic data. The impact of excluding 
this SN from this paper specifically and from the 3-year 
SNLS analyses more generally is negligible. 

The other two candidates were considered as "AGN?" 
from their real-time light curves and were not sent for spec- 
troscopy. The distances to their host galaxies are 0.2 and 0.3 
Tgai where Tgai is the effective galaxy radius introduced in 
Section These values are in the bulk of the distance dis- 
tribution, which peaks near 0.01 Tgai. There is no hint that 
these two events are highly centred, as would be expected 
for AGN's. Their light curves do not exhibit any peculiarity 
either when compared to normal SNe la. In both cases, the 
redshift of the associated host is in agreement with a su- 
pernova photometric redshift derived from the light-curves. 

Finally, SNLS detected eight supernovae which were 
classified as peculiar Type la after spectrosc o py, a s de- 
scribed inlBalland et alT (|2009t ). iBronder et all ()2008 f) and 
lEllis et all (120081^ : SNLS 03D3bb, SNLS 03D4cj, SNLS 
03D4ag, SNLS 03Dlcm, SNLS 04D3mk, SNLS 05Dlby, 
SNLS OSDlhk, SNLS 05D3gy. Su pernova SNLS 03D 3bb is 
a super-Chandrasekhar Type la ([Howell et al.l . 120061 ) with 
maximum light in the presurvey and does not pass the SN 
selection cuts. Supernova SNLS 03Dlcm is a 1991T-like 
object, with Xi — 4.54 and is therefore excluded by our 
cuts. Supernova SNLS 05Dlby is rejected because of its 
red colour (C = 0.66). The other five events do not exhibit 
any sign of peculiarity in their light curves and are in our 
photometric sample. 



4.3. Comparison with the SNLS 3-year cosmological sample 



4.4. Comparison with the SNLS sub-luminous SN la sample 



In lGuv et al.l (j2010[ ). SNLS describes a set of well sampled 
light curves for 252 Type la supernovae that was used to 
deriv e the 3-year cosmo logical analyses of the collabora- 
tion (jConlev et al.l . l201lt) . The photometric SN la selection 
described in this paper recovers 172 events of that sample. 
The lost fraction (32%) breaks down into 6% due to SN 
selections, 16% due to host galaxy redshift availability and 
10% due to SN la selections. Note that this result agrees 
with the efficiencies from synthetic SN la light curves that 
we quoted in Section B. 1.61 



iGonzalez-Gaitan et al.l (|201ll l provides a photometrically 
selected sample of 18 sub-luminous SNe la from SNLS data 
with z < 0.6. All but two of these SNe are included in the 
sample analysed in this paper. Among those 16 events, 7 
were kept by our photometric selection and 9 were rejected 
because of insufficient temporal sampling (4), cuts (2) 
or extreme values of SALT2 colour or Xi parameters (3). 
Despite the use of a light curve fitter trained on normal 
SNe la only, our photometric selection appears to have some 
efficiency on sub-luminous SNe la as well. 

4.5. Contamination by core-collapse supernovae 

Our set of 485 photometric SNe la candidates includes 175 
spectroscopically identified SNe la (see Table S]). One other 
event was spectroscopically identified as a possible core- 
collapse supernova initially, but was later found to be a 
SN la after it was discovered that the wrong spectrum had 
been extracted. The remaining 309 events have no conclu- 
sive spectroscopic typing either because no spectra were 
obtained or because of insufficient spectral signal-to-noise. 
In this section we will estimate how many of the 485 events 
are, in fact CC SNe that were incorrectly identified photo- 
metrically as SNe la. Core-collapse supernovae are typically 
more than a magnitude fainter than SNe la so a CC SN is 
unlikely to be confused with a SN la unless its redshift is 
significantly overestimated so that its faintness is explained 
by its distance. 

It is relatively simple to place an upper limit on the 
contamination at low-redshift, z < 0.4. Supernovae at these 
redshifts are spectroscopically unambiguous because of the 
visibility of the 615 nm Si II absorption feature characteris- 
tic of SNe la. While spectra for aU SNLS CC SNe were not 
obtained, those that appeared to be possible SNe la during 
the early phase of the explosion were given high priority 
for spectroscopy. We can therefore expect that for z < 0.4 
we have a complete spectroscopic sample of those CC SNe 
most likely to be photometrically selected as SNe la. Our 
photometric sample at z < 0.4 contains 35 spectroscopically 
identified SNe la and no spectroscopically identified CC SN, 
which gives a 95% CL upper bound of 3 core-collapse su- 
pernovae. The contamination of our sample at low redshift 
is thus expected to be less than ^ 3/35 ~ 9% at the 95% 
CL. This assumes that spectroscopic efficiencies at these 
redshifts are similar for Type la supernovae and for those 
among the CC events that would appear as possible SNe la. 

To refine this estimate and extend it to higher redshifts, 
we build a light curve model for CC SNe based on the SNLS 
sample of low-redshift {z < 0.4) supernovae, using the ana- 
lytical model of equation ^ to describe the light curves. To 
set up the simulation, we defined, as a CC test-sample, the 
set of 117 low redshift events used in lBazin et al.l (|2009[ ) to 
measure the CC rate at z < 0.4. The 117 events include the 
33 spectroscopically identified CC SNe and 84 events with 
no spectra and not photometrically selected as SNe la. To 
these, we added 9 spectroscopically identified CC SNe that 
had either no host galaxy photometric redshift or a host 
galaxy photometric redshift above 0.4. 

As detailed in section [31 the analytical model of equa- 
tion ID) describes the signal with four parameters in each 
filter, namely an amplitude, a date of maximum, a fall time 
and a rise time. Distributions of these parameters in the im 
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Fig. 16. Magnitude corresponding to Ai, the ZA^-band 
amplitude parameter of the model representing SN light 
curves, as a function of redshift, for fast-declining CC su- 
pernovae (top) and plateau ones (bottom), as measured 
in the 3-year SNLS data. Red dots come from fits to im- 
band light curves, green squares (resp. blue triangles) from 
fits in the tm (resp. (/m) filter taken as estimates of Ai at 
higher redshifts (see text). Symbols with black borders indi- 
cate spectroscopically identified events. The curves are fits 
to the average redshift dependence of the magnitude, ex- 
cluding points outside the dashed lines. As a comparison, 
black dots give magnitudes from spectroscopically identi- 
fied SNe la. 



filter, as well as correlations between these and the param- 
eters in each of the other filters were measured in the test- 
sample and parametrised for their subsequent use in the 
light curve simulation. This was done separately for events 
with plateau light curves (35 events, either events spec- 
troscopically confirmed as SNIIP or events with fall times 
above 60 days) and for fast-declining events (77 events, ei- 
ther events spectroscopically identified as SN lb, SN Ic or 
SN II or events with fall times below 40 days). No further 
distinction could be made in the latter category between 
Type Ib/c and Type II supernova light curves. 

The test-sample is mostly limited at low redshift, below 
z ~ 0.4. In order to model the average redshift depen- 
dence of the amplitudes, the parameters of the low-redshift 
supernovae were redshifted to mimic supernovae at higher 
redshifts. The amplitudes in the tm and gM bands at a low 
redshift z were thus taken as estimates of Ai at a higher 
redshift z,- defined as: 



l + Zi 



(1 + ^) 



after rescaling by: 
A, di(z)2 



where and A are the mean wavelengths of the and 
tm or gM filters and dL{z) is the luminosity distance for 
the cosmology introduced in Section 13 .3. II In the above for- 
mulas, z was taken as the spectroscopic redshift whenever 
available (50% of the cases), and as the host galaxy pho- 
tometric redshift otherwise. Amplitudes in the tm and gM 
filters were thus used to estimate the im band amplitude Ai 
at redshifts 0.4 < Zi < 0.7 and 0.7 < Zi < 1.2, respectively. 
Figure [11] presents the estimated Ai as a function of Zi . 

The average redshift dependence of the modelled u/- 
band magnitude was fitted with the logarithm of a third or- 
der polynomial in the redshift Zi. The test-sample of plateau 
(resp. other CC) supernovae has a scatter of 0.5 mag (resp. 
0.8 mag) around this average behaviour. Given the above 
scatters in magnitude, the dependence of our Ai modelling 
with the reference cosmology assumed in the luminosity dis- 
tance computation was found to be negligible. The fitted 
redshift dependence of the modelled ijvf-band magnitude 
and the spread of the residuals around the latter were used 
to simulate «j\/-band amplitudes according to a Gaussian 
distribution. Amplitudes in the other bands were simulated 
to reproduce the observed colours in the test-sample. The 
colours simulated at high redshift were checked to be identi- 
cal to those for lower redshift supernovae in the appropriate 
(bluer) bands. 




A,, = A- 



A dL{zt 



Fig. 17. Difference between redshift assignments as a func- 
tion of the SALT2 fitted colour, for synthetic core-collapse 
events (blue dots) and data events (all other symbols), af- 
ter all selections previous to that on SALT2 Xi and colour 
(see text). Green circles (resp. red triangles) are data events 
present in our sample which have been identified as Type 
la (resp. core-collapse) supernovae by spectroscopy. Open 
black circles stand for data events with no spectroscopic 
identification. 



We generated synthetic core-collapse light curves in the 
redshift range < z < 1.2. As in the SN la simulation, 
host galaxy photometric redshifts were also generated, in 
order to reproduce the selections of the analysis. A total 
of 40,000 light curves was simulated (20,000 in each cate- 
gory), of which 13% remained at the end of the SN selec- 
tion. To illustrate the output of the simulation. Figures [T7] 
and \TE\ show some of the constraints used in the SN la se- 
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Fig. 18. Colour-magnitude diagrams from SALT2 fitted 
magnitudes used in the photometric SN la selection, for 
synthetic core-collapse events and data events that have 
passed all selections prior to those based on these diagrams. 
Same colour code as in Figure [T71 



lection, with SNLS data now compared with synthetic CC 
light curves. We observe that our simulation based on low- 
redshift events reproduces the distribution of spectroscopi- 
cally identified CC SNe and of supernovae with no spectra. 

The contamination of the photometric sample of SN la 
candidates by CC SNe was then derived as follows. Our syn- 
thetic light curve simulation assumed redshift-independent 
volumetric supernova distributions. To account for this. 



synthetic events were first weighted by the actual red- 
shift dependence expected for the volumetric explosion 
rates, ass uming the SN la r ate to be proportional to 
(1 + z)2 (jPritchet et all . I1999D and the core-collapse SN 
rates to vary as the star for mation rate, i.e. as (1 -I- 
(|Hopkins fc Beacoml[2006[) . The numbers of weighted 
events passing the photometric SN la selections were then 
computed for the three samples of synthetic light curves 
(one for SNe la, one for SNe IIP and one for all other core- 
collapse SNe) in the redshift range z < 1.2. To normalise 
the SN CC and la simulations with each ot her, we used 
the CC to la volumetric rate ratio published in 'Bazi n et al.l 
(|2009D . that holds for low redshift {z < 0.4) events, m an 
absolute magnitude interval extending to 4.5 mags fainter 
than normal SNe la. The normalisation was set in order 
that weighted synthetic events reproduce the published ra- 
tio in the same conditions of redshift and magnitude. 

After normalisation, we obtained a contamination of 
17.6 ± 3.9{stat.) supernovae in our photometric sample of 
485 SN la candidates. Essentially all (99%) of the contam- 
inating events are non-plateau core-collapse supernovae. 
The statistical uncertainty is mostly due to the small num- 
ber of events in the sample of low-redshift CC SNe which 
led to the measurement of the CC to la volumetric rate 
ratio used in the above normalisation. The statistical un- 
certainty on that measurement was 20%. 



Table 5. Contamination from core-collapse supernovae. 



Redshift 


Events 


CC 


Zgal < 0.4 


49 


4.4± 1.1 


0.4 < Zgal < 0.8 


196 


11. 1± 2.5 


0.8 < Zgal 


240 


2.1± 0.5 



Notes. The table shows the number of events passing the pho- 
tometric SN la selection and the estimated CC contamination 
in three bins of host galaxy photometric redshift. 



Note that the synthetic CC events remaining at the end 
of the selections are at moderate host galaxy photometric 
redshift, Zgai ^ 0.5 on average with an r.m.s. of 0.2 (see 
also TableO and, when reconstructed as SNe la by SALT2, 
exhibit positive colours (with a mean value of 0.15 and an 
r.m.s. of the same order) and a broad distribution of Xi 
values (with a mean of -0.3 and an r.m.s. of 1.7). Moreover, 
the rate of 5a outliers in host galaxy photometric redshift 
is 5% in this sample, slightly higher than the 1.4% observed 
in the whole photometric sample. 

5. Comparison of event subsamples with and 
without spectroscopic identification 

The sample of 485 SN la photometric candidates contains 
175 spectroscopically identified SNe la, called hereafter the 
"identified sample", and 310 events with no spectral con- 
firmatior{3, called hereafter the "unidentified sample" . The 
distribution of the peak magnitude in im for the two sub- 
samples is given in Figure 1191 The unidentified sample is 



* Note that event SNLS 06D2bo has been included in the 
unidentified sample. 
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Fig. 19. Distributions of the SALT2 fitted iM peak mag- 
nitude (top) and of the host galaxy photometric redshift 
(bottom) for the identified (in green) and unidentified (in 
black) subsamples of photometrically selected SNe la. 



Table 6. Comparison of events with and without spectro- 
scopic identification. 



Parameter 


Identified 


Unidentified 


KS 




sample 


sample 


proba 


C 


0.00 ± 0.01 


0.02 ±0.02 


0.40 




0.23 ± 0.08 


-0.01 ±0.13 


0.30 



Notes. The mean values of SALT2 colour and Xi are given for 
photometrically selected events with jai < 23, split into subsam- 
ples with and without spectroscopic identification. Numbers in 
the last column are the Kolmogorov-Smirnov probabilities that 
the two distributions arise from the same parent distribution. 



The characteristics of the two subsamples of bright 
events were compared, based on their complete light curves 
as reconstructed in this analysis. Their colour and Xi dis- 
tributions were found to be compatible, as summarised in 
Table HI We also compared the colour-magnitude and Xi- 
magnitude relations in the two subsamples. For this pur- 
pose, the distance modulus of each event was defined from 
the SALT2 fitted B-band peak magnitude to|j, Xi and 
colour C as: 



M ± aXi - I3C 



(2) 



using values of A/, a and /3 introduced in Section 13.3.11 
Residuals were then computed from these distance mod- 
uli by subtracting 51og[P(2;)] where P{z) is a third-degree 
polynomial that we fitted on the total photometric sam- 
ple in order to describe the mean redshift dependence of 
the distance modulus with no assumption on a specific cos- 
mology model. The Hubble diagram residuals without the 
colour or Xi term in the distance modulus are represented 
in Fig. [201 The figure shows the fit of the colour- magnitude 
and Xi-magnitude relations from the full sample of bright 
events, as fits from the two subsamples were found to be 
indistinguishable. 

In conclusion, when restricted to the same range of im 
peak magnitudes, the unidentified and identified samples of 
events do not exhibit significant differences. 



about one magnitude deeper on average than the identi- 
fied sample, as expected from the constraints to obtain a 
spectrum for the identified SNe la. This translates into an 
average redshift of 0.60 for the identified set and of 0.86 for 
the additional set. 

5.1. Comparison of bright events 

At bright magnitudes, ia/ < 23, there are 110 identified 
events and 56 unidentified ones, all but 5 events found in 
the real time analysis. We traced back the reasons why 
these bright events missed a spectroscopic identification. A 
third of them had been detected in the real-time analysis, 
declared (from photometry) as probable SNe la and sent for 
spectroscopy. Insufficient quality of the spectrum, however, 
resulted in a spectroscopic redshift (usually from the host 
galaxy) but no typing of the supernova. Another third were 
also declared as probable SNe la but could not be scheduled 
for spectroscopic follow-up near maximum light. The last 
third were usually less convincing candidates according to 
their real-time, and thus partial, light curves. 



5.2. Comparison of full samples 

We now consider the full set of photometrically selected 
SNe la as an extension of the identified subsample towards 
fainter events, and estimate the impact on distance moduli 
of using the limited sample of spectroscopically identified 
SNe la. The Malmquist bias due to spectroscopic sample 
selections is an important issue in cosmology fits and was 
thoroughly studied in S NLS with Mon t e Car lo simulations. 
The results, reported in iPerrett et alJ (|2010t ) were used to 
correct SN la distance inoduli i n SNLS 3-year cosm olog- 
ical analyses (|Guv et all . 120101: IConlev et all l201lt) . The 
aim here is to check whether we can measure this bias di- 
rectly from data and with an analysis completely indepen- 
dent from that used to define the 3-year SNLS cosmological 
sample. 

To do so, ideally one would compare average distance 
moduli at a given redshift, measured from the identified 
subsample and from the whole photometric sample. Our 
statistical sample being limited, we compare Hubble dia- 
gram residuals in large bins of redshifts instead of distances 
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Fig. 20. Partial Hubble diagram residuals (and profiles) as 
a function of colour (top) and Xi (bottom) for the identi- 
fied (green filled circles and profiles) and unidentified (black 
open circles and profiles) subsamples of photometrically se- 
lected SNe la with zm < 23. The straight lines are linear 
fits to the full sample of events with < 23. 



at given redshift values. As all events were submitted to 
the same selection procedure and the same redshift assign- 
ment, such a comparison is directly sensitive to different 
Malmquist biases and should not be altered by systematic 
uncertainties related to data processing. 

The residuals of the Hubble diagram are presented in 
Fig. [5T] for the two samples. The residuals were again 
defined from the distance moduli in Eq. [5] by subtract- 
ing 51og[P(z)] where P{z) is the third-order polynomial 
introduced in the previous section. The Hubble diagram 
residuals were then averaged in three large redshift bins, 
[0.3 - 0.6], [0.6 - 0.8] and [0.8 - 1.05], where both identi- 
fied and unidentified samples had enough events to allow a 
quantitative comparison. In the computation of the means, 
3(7 outliers were rejected. 

Mean values of the residuals (each weighted by the in- 
verse of its variance) are reported in Table [T] The vari- 
ance on each point includes both a contribution from the 
light curve fit and one from the use of a photometric red- 
shift. For the latter, we used a typical redshift uncertainty 
of Sz — 0.05, that we translated into a distance modulus 



Fig. 21. Hubble diagram residuals. Green dots stand for 
the identified SN la subsample and open black circles for 
the unidentified subsample. Means are shown in three red- 
shift bins, black for identified events, red for full sample. 
The red upward triangle is event SNLS 06D2bo which was 
erroneously typed as a possible core-collapse supernova ini- 
tially (see Section IT2l for details). The pink crosses mark 
the two 3(7 outliers, which were excluded in the computa- 
tion of the means. 



uncertainty using the derivative of the cosmological model 
introduced in Section 13.3.11 By construction, the residuals 
of the total sample are (see also Figure I^Tj) . The last line 
of Table [7] gives the difference between the mean residuals 
in the identified and total samples, which directly measures 
the mean difference in distance modulus between the two 
samples, Afiid- 

Up to a redshift of 0.8, distance moduli in the identified 
and the total samples do not exhibit significant differences. 
In the last redshift bin, corresponding to a mean redshift 
of 0.9, distance moduli in the two samples have an average 
offset of Afi,d = -0.030±0.036 (stat.). This resuh is stable 
with respect to different redshift binning or outlier cuts. 

5.3. Offset correction and systematics 

As shown in Figure [22j the bias towards selecting bluer, 
and hence brighter events at high redshift is clear in the 
identified sample and attenuated in the photometric sam- 
ple. Although the latter effect may be due to evolution of 
Type la supernovae, it can also result from selection effects 
(our photometric sample being magnitude limited and bi- 
ased against SNe in faint host galaxies). In that case, the 
offset value quoted in the previous section should be cor- 
rected for the incompleteness of the photometric sample. A 
correction for the residual core-collapse contamination of 
the sample may also be needed. 

These two corrections were determined from the syn- 
thetic supernova light curves, which were submitted to 
the same selections as data. Residuals were here com- 
puted w.r.t. the luminosity distance for the cosmology 
used in the simulation (see Section I3.3.ip . In the red- 
shift bin [0.80, 1.05], synthetic SNe la passing our pho- 
tometric selection exhibit a negative average residual of 
-0.0148 ± 0.006(stat.), while the selected CC SNe (less 



19 



G. Bazin et al.: Photometric selection of Type la supernovae in the SNLS 
Table 7. Hubble diagram mean residuals for different samples of SN la data selected by this analysis. 



SN la sample Mean residuals in redshift bins 

0.3 < z < 0.6 0.6 < z < 0.8 0.8 < z < 1.05 

Identified 0.004 ± 0.032 (76) 0.003 ± 0.033 (51) -0.025 ± 0.041 (31) 

Unidentified -0.002 ± 0.053 (31) -0.009 ± 0.032 (67) 0.013 ± 0.022 (160) 

Total 0.002 ± 0.027 (107) -0.003 ± 0.023 (118) 0.005 ± 0.019 (191) 

Anid = Identified - Total 0.002 ± 0.017 0.006 ± 0.024 -0.030 ± 0.036 

Notes. Host photometric redshifts were used for all events. The number of events in each bin is indicated in parenthesis. The last 
line shows the difference between the mean residuals in the identified and total samples. Uncertainties are statistical only. 
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Table 8. Effect of redshift and photometry resolutions on 
mean residuals. 



Analysis 
conditions 


Mean residual 
Bright SNe Photometric SNe 


This analysis 


-0.033 ± 0.007 


-0.015 ±0.006 


\Az\/(l + z) < 10% 
\Az\/{l + z) <3% 
\Az\/{l + z) < 1% 

Perfect photom. & z 


-0.034 ± 0.007 
-0.042 ± 0.009 
-0.044 ± 0.014 
-0.043 ± 0.007 


-0.016 ±0.006 
-0.020 ± 0.008 
-0.021 ±0.012 
-0.024 ± 0.006 


Difference 


-0.010 ± 0.009 


-0.009 ± 0.008 



Notes. The table details the changes in the weighted mean 
residuals of the identified and photometric SN la samples, as ex- 
pected if resolutions in redshift and photometry were gradually 
improved w.r.t. their values in the present analysis (first line). 
All changes were estimated from synthetic SN la light curves. 
Uncertainties are statistical. 
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Fig. 22. Distributions of colour (top) and Xi (bottom) as 
a function of the host galaxy photometric redshift, for the 
identified (green filled circles and profiles) and unidentified 
(black open circles and profiles) subsamples of photometri- 
cally selected SNe la. 



than 1 event out of 191) contribute for a tiny residual of 
+0.0005 ± 0.0008(stai.). The offset of the photometric 
sample thus amounts to Afiphot = —0.014 ± O.OOQ{stat.) 
and the total offset of the identified subsample in the red- 
shift interval 0.8 < z < 1.05 is A/i = A/i^d + Ajiphot = 
-0.044 ±0.036(stoi.). 

Systematic uncertainties due to the limited precision 
of redshift assignment and photometry may affect this re- 



sult. These effects were again studied with the synthetic 
SN la light curves. Changes in the mean residuals of the 
photometric sample and its bright subsample are reported 
in Table [HI The difference between the residuals of the two 
samples compares directly with A/i^d, while the mean resid- 
ual of the photometric sample identifies with the correction 
term A^phot- Table [5] shows that improving the redshift de- 
termination and the photometry precision leads to almost 
identical shifts of the mean residuals of the two samples. 
These two effects have therefore a negligible impact on the 
measurement of Afiid- The table also shows that the sys- 
tematic uncertainty on A^phot is of order -0.01, negligible 
w.r.t. the statistical uncertainty of the method. 

As previously mentioned, a precise determination of the 
magnitude offs et in the SNLS 3-yea r spectroscopic sample 
was obtained in' Perrett et al.l (|2010t l with Monte Carlo sim- 
ulations. The shift in the average magnitude of the spectro- 
scopic sample towards brighter values was found to become 
significant for z > 0.75, rising from — 0.013 ± 0.001 (stat.) 
at z ~ 0.8 to -0.038 ± 0.003(stat.) at z - 1.05, with a 
systematic uncertainty of 20%. The result reported here, 
Afj. = -0.044 ± 0m6{stat.) ± 0.010(syst.) in the redshift 
interval 0.8 < z < 1.05, obtained mostly from data and 
from an analysis independent from the real-time SNLS pro- 
cessing, agrees with the above values. 



6. Conclusions 

The SuperNova Legacy Survey offers one of the best data 
sets to test SN la photometric typing methods in the red- 
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shift range between 0.2 and 1.0. Although spcictra remain 
essential to achieve the level of accuracy required in cosmol- 
ogy studies, many supernova studies such as rate measure- 
ments or correlations with host galaxy properties benefit 
from the larger statistics allowed by photometrically iden- 
tified samples. 

In this paper, the 3-year SNLS image sample was sub- 
mitted to a deferred search for Type la supernovae based 
on their multi-colour light curves and an external catalogue 
of photometric redshifts of host galaxies. SN-like transient 
events were first selected rising light curve shape criteria. 
The selected light curves were fitted under the assumption 
that the events were SNe la using SALT2 as light curve 
fitter. Selections were then designed to both discriminate 
SNe la against core-collapse ones and to assess the relia- 
bility of the redshift assignment. They were applied to the 
SALT2 fitted parameters, rest-frame B — V colour, stretch- 
related Xi parameter, as well as colours and magnitudes in 
the MegaCam filters. These selections allowed us to select 
a sample of 485 SN la events with host galaxy photometric 
redshifts. 

The performance of the analysis was studied with syn- 
thetic light curves of Type la and core-collapse supernovae. 
The selection efficiency for bright SNe la with well sampled 
light curves is 80% in case a host galaxy photometric red- 
shift is available. The contamination of the selected sample 
from other types of supernovae is 4%. 

One third of the photometrically selected events have 
confirmed la types from SNLS spectroscopy. These events 
were compared with the rest of the photometric sample 
which has no or incomplete spectroscopic information. In 
the range of magnitudes allowed by spectroscopic identi- 
fication, there is no significant differences in the intrinsic 
features of events with and without spectroscopic identifi- 
cation. 

The higher magnitude limit of the whole photomet- 
ric sample allowed distances derived from spectroscopi- 
cally identified events to be tested in order to measure the 
Malmquist bias due to spectroscopic sample selections, di- 
rectly from data and with an analysis independent from the 
SNLS cosmological analyses. A magnitude offset was found 
at redshifts above 0.8, though with a large statistical un- 
certainty. The obtained offset is consistent with the more 
precise result derived from Monte Carlo simulations which 
was used to correct SN la distance moduli in the SNLS 
3-year cosmological analyses. 

This paper demonstrates the feasibility of a photometric 
selection of high redshift supernovae with known photomet- 
ric host galaxy redshifts, opening interesting prospects for 
cosmological analyses from large photometric SN la sur- 
veys, as planned in the near future. In this respect, the 
points we found most important for such an analysis are 
the following. The cadence of observations must be chosen 
to obtain light curves with a temporal sampling of a few 
days. The number of filters must allow SN colours and mag- 
nitudes to be determined in a reliable way on the whole SN 
redshift range. A low contamination by core-collapse super- 
novae can be achieved when using photometric host galaxy 
redshifts with a few % resolution and rate of catastrophic 
redshifts, provided selections include enough criteria to re- 
ject events with incorrect redshift assignments. 
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Appendix A: detection efFiciency model 

The detection efficiency for a supernova with date ipi of 
maximum Hght in iM-band depends on its peak magnitude 
moi, on the seeing and sky background during nearby ob- 
serving times tk, and on the relative epochs toi — t^. The 
dependence of the efficiency on these variables was studied 
with Monte Carlo images on Dl. 
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Fig. 23. Detection efficiency e(moi) on Dl derived from 
Monte-Carlo images as a function of the generated magni- 
tude moi at maximum light. Blue diamonds illustrate the 
efficiency for the first season of data (averaged over the full 
mosaic excluding CCD 3), red disks are for the second sea- 
son (full mosaic), open black circles are for the third season 
(full mosaic) , and pink triangles are for the supernovae oc- 
curring during the first season on CCD 3. 



As an example, the season-averaged average detection 
efficiency e(moi) is illustrated in Figure [531 On average over 
the CCD mosaic, all three seasons exhibit very similar effi- 
ciencies, except for CCD 3 which failed to work during three 
out of the six dark time periods of the first season. The effi- 
ciency is nearly magnitude-independent out to moi = 23.5 
and then steeply declines at fainter magnitudes. 

With varying observing conditions, we found that the 
efficiency as a function of magnitude had the form of the 
mean efficiency in Figure 1231 but that the maximum effi- 
ciency and cutoff magnitude depended on conditions. The 
adopted efficiency model therefore has the form: 

e(toi) = f-max{hi) F{m - rndtoi)) 



— 0.9 




0.6 0.8 1 1.2 1.4 

seeing - 0.55 (background - 20) 



Fig. 24. Lunation detection probabilities a{oi) from the 
Dl Monte Carlo as a function of mean lunation seeing (in 
arcsec) and background (in mag per arcsec^). The line is 
the adopted a{oi). 
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Fig. 25. Analytical model of the temporal dependence 
f3{Ati), here illustrated for a fourteen-day lunation. Dashed 
lines indicate the lunation boundaries. The model consists 
of two decreasing exponentials outside the lunation bound- 
aries (with a steeper decline for positive Ati because the fall 
time of a SN la is larger than its rise time) and is roughly 
flat inside. 



where F{m — rric) is the mean of the functions shown in 
Figure [12] (F(-mc) = 1 and F(0) = 1/2). The function 
£max is the efficiency at bright magnitudes while rric is the 
magnitude at which the efficiency drops to tmaxl^- The t^i 
dependence of tmax and vfic reflects the observing condi- 
tions near ioi- It was found that they could both be deter- 
mined by a single function /(toi)" 

emax(^o^) = 1-0.2/(^0,:) m^{t^,) = 25.3-5.4/(to0 (3) 

The efficiency- loss function, f{tQi), vanishes for the best 
observing conditions. It was taken to be a product over all 
lunations: 

^ll[l-aioi)x(3{AU)] (4) 
I 

where a is a function of the mean observing conditions, 
oi = [seeing, sky background) for the lunation and /3 is 
a function of the difference, Ati, between the time of su- 
pernova maximum light and the time of the nearest mea- 
surement. Figures [m and [55] show the adopted functions. 
Note that for Ati > 50 days, /3(At;) — > 0, so in practice 
only the three lunations closest to the explosion determine 
f{toi)- The impact of temporary CCD failures is automat- 
ically accounted for in this model by setting a{oi) — for 
non-working CCDs, leading to large intervals between the 
maximum light date toi and the actual dates of observation, 
and hence to lower detection efficiencies. 



22 



G. Bazin et al.: Photometric selection of Type la supernovae in the SNLS 



o r I X 1 r 

E - .• • i : . 

25 - u - •rfe»3f 

0.98 - • •V'°o' S 

24.75 - o o 

■ » o ■ 0.96 - * S 1 

on 
24.5 7 ■ 

0.94 - 

24.25 - 

o 

0.92 - 

24 - 

- ' 0.9 r 

0.1 0.2 0.1 0.2 0.3 

f(»Oi) *(toi) 

Fig. 26. Correlations between the measured parameters 
(rric and Emax) and /(toi)i as observed in the Dl Monte- 
Carlo. Each point is an average over all toi dates in a lu- 
nation (filled symbols) or inter-lunation period (open sym- 
bols). Circles are averaged over the full mosaic, except for 
the first season where CCD 3 has been singled out (red 
squares) due to intermittent failures. Linear fits are used in 
the model to describe the observed correlations. 



to show all cases of inter-lunation durations, from fifteen 
days (normal CCD working conditions) to over one month. 
Figure [25] shows that the observed correlations can be de- 
scribed by the same linear fits, whether ioi occurred during 
a lunation or outside, for working CCDs or not. 

Figure E71 compares the model for 1 — /(toi) (blue curve) 
with the 1 — f{toi) derived directly from the Monte-Carlo 
images (shown as the black points). There is a clear agree- 
ment between the two estimates. The model reproduces 
both the efficiency during the observation periods and the 
loss of efficiency between these, which varies with the du- 
ration of the absence of data. This is clearly emphasized in 
the bottom plot of Figure [57] where the temporary failures 
of CCD 3 enhance the efficiency loss due to lack of data. 

Appendix B: additional efFiciency plots 

This appendix gives the SN la photometric selection effi- 
ciency as a function of the generated SN la colour and Xi . 
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Fig. 27. Temporal dependence of the detection efficiency, 
l~ f{toi), on Dl, with dates originating at January 0, 2003. 
The model (blue curve) is compared with Monte-Carlo es- 
timates (black points) . Actual observation dates are shown 
as the lower red histograms. The bottom plot is for CCD 
3 which failed to work during the second and the last two 
lunations. The top plot is for the rest of the mosaic. 
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The linear dependence of emax and nic on f{tQi) (see 
equation [31) is illustrated in Figure [21] We see that the pri- 
mary effect of degraded observing conditions is to lower 
the cutoff magnitude: conditions that degrade rric ~ 25 to 
TUc 24 change the efficiency for bright supernova only 
from emax ~ 0.98 to e^nax ^ 0.96. Note also that nic differs 
on average by about 0.25 mag for supernovae occurring dur- 
ing a lunation compared to those between lunations. The 
case of CCD 3 in the first season has been also singled out 



Fig. 28. Selection efficiency from synthetic SN la light 
curves as a function of the generated colour and Xi, at 
different stages of the analysis. 
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